Ab-initio study of the vibrational properties of Mg(AlH4)2 
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Based on Density Functional Theory and Density Functional Perturbation Theory we have stud- 
ied the thermodynamical and vibrational properties of Mg(AlH4)2. The crystal structure recently 
proposed on the basis of x-ray powder diffraction data has been confirmed theoretically by the com- 
parison of the experimental and theoretical IR and Raman spectra. The main discrepancy regards 
the position of the hydrogen atoms which makes the theoretical AIH4 tetrahedra more symmetric 
than the experimental ones. The calculated thermodynamical decomposition temperature is also in 
good agreement with experimental result. 



I. INTRODUCTION 

In the search for suitable materials for reversible hy- 
drogen storage the class of metal alanates has attracted 
much attention in recent years due to their high hy- 
drogen content, but primarly because of the possibil- 
ity to accelerate the kinetics of hydrogen uptake by 
doping^'^. Sodium alanate, the prototypical material of 
this class, has long been known. Although its thermo- 
dynamics decomposition temperature with hydrogen re- 
lease is relatively low, the kinetics of the dehydrogena- 
tion/rehydrogenation reactions is very slow. The interest 
in this material for hydrogen storage increased dramati- 
cally in 1997 when Bogdanovic et aii demonstrated that 
reversible storage can be achieved at moderate temper- 
ature and hydrogen partial pressure by the addition of 
catalysts, most notably titanium. This discovery boosted 
an intense activity on the study of other members of the 
alkali metal alanates as well^'- . 

More recently, magnesium alanate Mg(AlH4)2, as a 
representative of the alkali earth metal alanates, has 
also been considered as a possible materials for hydro- 
gen storage^iSi^. Mg(AlH4)2 decomposes readly in the 
temperature range 110 °C - 200 °C according to the re- 
action 

Mg{AlHi)2 MgH2 + 2AI + 3F2 (1) 

which would correspond to a maximum reversible hy- 
drogen content of 7 wt%. Although the reversibility of 
reactionQlhas not been demonstrated yet, the easy of the 
decomposition reaction even in the absence of doping, 
suggests that, in analogy with NaAlH4, the insertion of 
suitable catalyst might enhance the kinetics of hydrogen 
uptake. The structure of Mg(AlH4)2 has been recently 
resolved from X-ray powder diffraction pattern aided by 
quantum-chemical calculations on cluster models which 
allowed to tentatively assign also the position of hydrogen 
atoms, not detectable by X-ray diffraction*. 

In this work, we present an ab-initio study of the ther- 
modynamical and vibrational properties of Mg(AlIl4)2 
aiming at providing a better estimate of the decompo- 
sition temperature from calculations on periodic models 
and at confirming the structure inferred experimentally 



from the comparison of theoretical and experimental IR 
and Raman spectra. 



II. COMPUTATIONAL DETAILS 

Calculations are performed within the framework of 
Density Functional Theory (DFT) with a gradient cor- 
rected exchange and correlation energy functional^, as 
implemented in the codes PWSCF and PHONONSiS. 
Norm conserving pseudopotentialaii and plane wave ex- 
pansion of Kohn-Sham (KS) orbitals up to a kinetic cutoff 
of 40 Ry have been used. Non linear core corrections are 
included in the pseudopotential of magnesiumi^. Bril- 
louin Zone (BZ) intergration has been performed over 
Monkhorst-Pack (MP)!^ 6x6x6, 4x4x4 and 16x16x16 
meshes for Mg(AlH4)2, MgH2 and metallic Al, respec- 
tively. Hermite-Gaussian smearing^"^^. of order one with a 
linewidth of 0.01 Ry has been used in the reference calcu- 
lations on metallic Al. Equilibrium geometries have been 
obtained by optimizing the internal and lattice structural 
parameters at several volumes and fitting the energy ver- 
sus volume data with a Murnaghan functior^iii. Residual 
anisotropy in the stress tensor at the optimized lattice 
parameters at each volume is below 0.6 kbar. Infrared 
and Raman spectra are obtained from effective charges, 
dielectric susceptibilities and phonons at the F point 
within density functional perturbation theoryi^. Since 
Mg(AlIl4)2 is an uniaxial crystal, the infrared absorption 
spectrum depends on the polarization of the trasmitted 
electromagnetic wave with respect to the optical axis. 
Relevant formula for the calculation of the IR and Ra- 
man spectra for a polycristalline sample to be compared 
with experimental data are given in section II. 



III. RESULTS 
A. Structure and energetics 

The structure of Mg(AlIl4)2 has been recently assigned 
to a Cdl2-like layered crystal from X-ray powder diffrac- 
tion datai^. Rietveld refinement assigned the space group 
Piml (-D|rf) and lattice parameter a=5=5.199 A and 
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TABLE I; Experimental and theoretical (in parenthesis) po- 
sitions (in crystal units) of the four independent atoms at 
the experimental equilibrium lattice parameters (space group 
P3ml, a=5.199 A, c=5.858 A). The Wyckoff notation is used. 
If not reported the theoretical positions coincide with the ex- 
perimental ones by symmetry. HI and H2 are non-bridging 
and bridging hydrogen atoms, respectively. 



Mg(la) 











Al(2d) 


1/3 


2/3 


0.7 (0.7053) 


Hl(2d) 


1/3 


2/3 


0.45 (0.4349) 


H2(6i) 


0.16 (0.1678) 


-0.16 (-0.1678) 


0.81 (0.8111) 



c=5.858 A. The position of the four independent atoms 
in the unit cell are given in Table I. In the AIH4 tetra- 
hedron there are two Al-H bond lengths: the length of 
the Al-Hl bond with non-bridging HI hydrogen aligned 
with the c axis and the length of A1-H2 bond with the 
H2 atom bridging between Al and Mg. 

A picture of the crystal structure is given in Fig. ^ It 
can be seen as a stacking along the c axis of A1H4-Mg- 
AIH4 neutral trilayers. 




FIG. 1: (Color online) Trigonal structure of Mg(AlH4)2. The 
unit cell containing a single formula unit is drawn. 

The theoretical equilibrium lattice parameters ob- 
tained from the Murnaghan equation of state are a=5.229 
A (exp. 5.199) and c= 6.238 A (exp. 5.858). In the ge- 
ometry optimization we have also relaxed the constrained 



of P3ml symmetry, but we have always recovered the 
experimental space group. The misfit in the c axis ( 6 
%) is larger than usual in DFT-based calculations. As 
already pointed out in previous cluster calculations on 
Mg(AlH4)2^ the trilayers stacked along the c axis are 
neutral and not chemically bonded each other. Thus 
van der Waals interactions have been proposed to play 
an important role in interlayer cohesion. Since van der 
Waals interations are not present in current approxima- 
tions to the exchange and correlation energy functional, 
the interlayer bonding is underestimated which implies 
an expansion of the c axis with respect to the experi- 
mental data. The energy difference between the relaxed 
configurations at the experimental and theoretical equi- 
librium lattice parameters is nevertheless small ( 38 meV 
per formula unit). The theoretical atomic positions at 
the experimental lattice parameters are given in Table 
I. The theoretical geometry of the AIH4 tetrahedron is 
closer to an ideal tetrahedron than the experimental ge- 
ometry. In fact the theoretical length of the two Al-H 
bonds are A1-H1= 1.584 A (exp. 1.46) and Al-H2= 1.614 
A (exp. 1.69). The Al-H bond lengths do not change siz- 
ably (within 0.01 A) by changing the lattice parameters 
from the experimental values to the theoretical equilib- 
rium values. 

The electronic band structure along the high symmetry 
directions of the Irreducible Brillouin Zone is reported in 
Fig. 12 The analysis of the projection of the KS states 
on the atomic orbitals at the F point shows that the 
top of the valence band is mainly formed by 3p states 
of aluminum and Is state of hydrogen while the bottom 
of the conduction band is mainly formed by 3s states 
of magnesium. The band gap is thus a charge transfer 
excitation. 




A r MKT 



FIG. 2: Electronic band structure along the high symmetry 
direction of the Brillouin Zone. High symmetry points are 
labeled following Zak (Reff^). The zero of energy is the top 
of the valence bands. 

By neglecting at first the vibrational contribution to 
the entropy of the solids, the thermodynamical decom- 
position temperature Td of reaction ^ can be estimated 
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Td 



AE + PAV 
3(5'h2 — R) 



(2) 



where is the difference in volume of the sohd reac- 
tants and products of reaction ^ P is the pressure, AE 
is the energy difference of reactants and products in ^ 
S'hj is the entropy per mole of gaseus H2 and the gas 
constant R in the denominator comes from the PV term 
of gaseus H2. AV and AE are obtained from the ab- 
initio calculations. For MgH2 the theoretical equilibrium 
structural parameters are a— 4.470 A (exp. 4.50li^), c= 
2.943 A (exp. 3.010Ji) and the 4/ position (in Wyckoff 
notation of the P42/mnm space group) of the indepen- 
dent H atom in the unit cell is (0.3046,0.3046,0) (exp. 
(0.3040,0.3040,0)1^). The theoretical equihbrium lattice 
parameter of metalhc Al is a= 4.065 A (exp. 4.050 A2i). 
From the total energies at equilibrium we obtain AE— 
144 kJ/mol (48 kJ/mol per H2 molecule) which becomes 
AE= 139.5 kJ/mol (46.5 kJ/mol per H2 molecule) by 
including the zero point energy of the optical phonons 
(at the r point only) and of the acoustic bands within 
a Debye modeli^. This result is in good agreement with 
previous estimate from the extrapolation of cluster calcu- 
lations (123 kJ/mol^, 41 kJ/mol per H2 molecule). The 
PAV in Eq. |21 is negligible (4.4 J/mol). By assuming 
S'h2= 130 J/mol at 300 KiS. and atmospheric pres- 
sure, we obtain T]j= 111 °C. The predicted T^i shifts 
by 2 °C by including the temperature dependence of Sn^ 
as given by the Sackur- Tetrode expression^*^ , the vibra- 
tional contribution to the free energy of the solids and the 
translational and rotational energy of gaseous H^l^. The 
theoretical Td is comparable with the experimental de- 
composition temperature which falls in the range 110 °C 
- 200 °C which implies that kinetics effect do not affect 
dramatically the decomposition reaction of Mg(AlH4)2. 



B. Vibrational properties 

Phonons at the F-point can be classified according to 
the irreducible representations of the D^d point group of 
Mg(AlH4)2 as F= A Aig + 2 + 5Eg + 5 A2„ + 6 E„. 
Aig and Eg modes are Raman active while A211 and E„ 
modes are IR active. One and one E„ mode are uni- 
form translational modes. The IR active modes display 
a dipole moment which couple to the inner macroscopic 
longitudinal field which shifts the LO phonon frequen- 
cies via the non-analytic contribution to the dynamical 
matrijii^ 



D^^(n n') ^ '^^ ^c,,a'{K)qa'Zp^f3,(K')qp, 



Vo 



q • £^ 



where Z and e°° are the effective charges and electronic 
dielectric tensors, Vo is the unit cell volume and q is the 



phononic wavevector. The effective charge tensor for the 
three independent atoms in the unit cell (cfr. Table P) 
are 



/ 2.180 



ZiMg) = 



2.180 



(4) 



1.974 



1.934 

Z{Al) = I 1.934 



(5) 



1.678 



-0.627 
Z{m) = I -0.627 



(6) 



-0.641 



/ -0.564 
Z{H2) = -1.029 0.248 

\ 0.289 -0.672 y 

The electronic dielectric tensor is 



(7) 



2.894 



2.894 



(8) 



2.781 



The calculated phonon frequencies at the F point, ne- 
glecting the contribution of the longitudinal macroscopic 
field (Eq. ^ , are reported in Table II for the equilibrium 
geometry at the experimental lattice parameters. The 
phonon frequencies at the theoretical lattice parameters 
differ at most by 10 cm^^ with respect to those reported 
in Table II. 

For a uniaxial crystal like Mg(AlIl4)2, the macroscopic 
field contribution to the dynamical matrix (Eq. ^ in- 
troduces an angular dispersion of the phonons at the F 
point, i.e. the limit of the phononic bands cj(q) for q — s- 
depends on the angle 9 formed by q with the optical axis. 
The angular dispersion of the A2U and E„ modes due to 
the macroscopic field is reported in Fig. O 



1. IR spectrum 

Tha absorption coefficient of an uniaxial crystal de- 
pends on the the polarization of the trasmitted light with 
respect to the optical axis. The optical properties of the 
crystal can be obtained from the two dielectric functions 
e±{oj) and e|| {lu) which represent the response of the crys- 
tal to electromagnetic wave with electric field perpendic- 
ular (E _L c) and parallel (E || c) to the optical axis, 
respectively. The trasmitted electromagnetic wave at a 
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TABLE II: Theoretical phonon frequencies at the F point, os- 
cillator strengths (/j in Eq. IIUII of IR it- modes and coefficients 
of the Raman tensor of the Raman active mode, a, b for Aig 
and c, d for Eg modes (in unit of 10"^ Vo= 0.2743 A^, see 
section III B.2). The contribution of the inner longitudinal 
macroscopic field is not included (LO-TO splitting) (see Fig. 
^for the displacement patterns). 



Modes 


Energy (cm ^) 


fj 


a (c) 


b (d) 


Eff (1) 


87 




0.976 


0.578 


A29 (1) 


169 








A19 (1) 


232 




2.337 


1.210 


E„ (1) 


282 


3.162 






Eff (2) 


298 








A2u (1) 


302 


0.103 






A29 (2) 


355 








E„ (2) 


360 


0.758 






Eu (3) 


620 


3.896 






A2„ (2) 


663 


0.784 






E„ (4) 


716 


0.064 






E9 (3) 


742 




-0.457 


3.653 


E9 (4) 


758 




2.746 


0.023 


A19 (2) 


812 




3.970 


0.681 


A19 (3) 


1845 




8.768 


17.540 


A2„ (3) 


1850 


0.072 






Eg (5) 


1852 




5.202 


0.339 


E„ (5) 


1905 


0.597 






A2u (4) 


2013 


0.040 






Ai. (4) 


2077 




12.270 


5.160 
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FIG. 3: Angular dispersion of the F-point phonons. Only A2U 
and Eu modes have angular dispersion. 9 is the angle formed 
by the phonon wavevector and the optical axis. 



dinary wave with electric field lying in the plane formed 
by the optical axis and q. The dielectric function e±(uj), 
independent on 9, describes the response to the ordinary 
waves while the dielectric function ee((x') for the extraor- 
dinary wave is ^-dependent and is given by 



ee (uj) 



e_L(w)e||(w) 



e±(Lu)sin'^9 + e|| {uj)cos^ 



(9) 



The extraordinary wave coincides with the ordinary 
wave for 6 — 0, i.e. propagation along the optical axis. 

i±{uj) and e||(a;) can be obtained from ab-initio 
phonons, effective charges and electronic dielectric ten- 
sor as 




N 




(11) 



where u and ^ are the number of E„ and TO 
modes, respectively. The phonons entering in Eq. [TUl lfTT)) 
have wavevector q parallel (perpendicular) to the optical 
axis. The notation || and _L refers to the orientation with 
respect to the optical axis of the dipole moment of the 
phonon which coincides with that of the electric field of 
the trasmitted wave it couples to. The sum over k run 
over the N atoms in the unit cell with mass M^. e{j, k) 
and LOj are the eigenstates and eigenvalues of the dynam- 
ical matrix at the T point, without the contribution of 
the macroscopic field which has no effect on the purely 
TO modes. The absorption coefficient for the ordinary 
wave is given by 



— Imej^(a; + 17, 7 ■ 
nc 




]=1 K=l 



(12) 



c is the velocity of light in vacuum and n is the refrac- 
tive index. The absorption coefficient for the extraordi- 
nary waves Q!||(w) with E || c (and q _L c) is given by the 
analougous expression by changing v with /i in the sum 
over phonons in Eq. 1121 The (5-functions in Eq. 12 are 
approximated by Lorenztian functions as 



generic wavevector q forming an angle 9 with the op- 
tical axis would split in an ordinary wave with electric 
field perpendicular to the optical axis and in an extraor- 



S{UJ - UJj) = - 



2 



(13) 
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FIG. 4: (Color online) Displacement pattern of the IR and 
Raman active intra-tetrahedra modes. The lattice modes at 
lower frequencies are not reported. 



The functions a±{tu) and a||(a;) are shown in Fig. 
The phonons around 700-800 cm~^ are bending modes 
of the Al-H bonds. The highest frequency modes in the 
range 1800-2000 cm^-'^ are stretching modes of the Al-H 
bonds. In particular the A2ti(3) mode is a stretching of 
the Al-Hi apical bond. The other modes are stretching 
of the AI-H2 bridging bonds. The displacement pattern 
of the IR and Raman active modes are shown in Fig. 0] 
The modes below 360 cm~^ (cfr. Table ITl)l are lattice 
modes involving a rigid motion of the teatrahedra. 
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FIG. 5: IR absoption spectrum for a) electric field E || c (a||) 
and for b) E _L c (a±). The Lorenztian broadening is 7= 15 
cm-i (Eq. [SJ. 



Since the experimental data are available only for a 
polycrystalline samplel an angle-averaged absorption co- 
efficient is needed to compared theoretical and experi- 
mental data. For a generic q forming an angle 9 with the 
optical axis the absorption coefficient for the extraordi- 
nary wave can be obtained from the imaginary part of 
ee (Eq. O. For non-polarized light the total absorption 
coeffient can be equivalently expressed as 

"«H = lbEiE^Ai-^m--^.), (14) 

where e(j, k) and Ooj are eigenstates and eigenvalues 
of the full dynamical matrix including the non-analytic 
term (Eq. O which mixes E„ and A2M modes. 

The angle-averaged absorption coefficient has been ob- 
tained from Eq. 1141 as 

aave = ^ sin{9n)ae„ , (15) 

n 

where the sum runs over ten angles equally spaced 
in the range O-tt. The resulting theoretical absorption 
coefficient for a polycrystalline sample is compared with 
the experimental IR spectrum in Fig. |S1 




wave number (cm ') 

FIG. 6: (Color online) IR absorption spectrum for polycrys- 
talline Mg(AlH4)2. a) Theoretical spectrum with Lorentzian 
broadening 7= 15 cm~^ (Eq. 1131 . b) Theoretical spectrum 
(black line) with 7= 40 cm~^ compared with the experimental 
spectrum (gray line (red online), adapted from Ref.»). The ex- 
perimental spectrum is normalized by setting equal the exper- 
imental and theoretical intensity of the highest peak around 
1800 cm-\ 

Good agreement with experimental data is obtained 
with a Lorentzian broadening of 7= 40 cm~^. (Fig. 
EJl, but for the experimental shoulder around 800 cm~^ 
which is absent in the theoretical spectra. This misffi 
might be either due to an underestimation of the A2U (2) 
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mode which is very close to the strongest IR mode E2„(3) 
(cf. Fig. O or to a large inhomogeneous broadening of 
the experimental spectra. In fact, the line- width of the 
experimental peak are very large and might be partially 
due to residual solvent adducts which are released only 
after dehydrogenationi. 

2. Raman spectrum 

The differential cross section for Raman scattering 
(Stokes) in non-resonant conditions is given by the fol- 
lowing expressionSiiSS (for a unit volume of scattering 
sample) 

= E ^l^'s • • eLp(n(c.) + l)5{u: - c,), (16) 
j 



The experimental Raman spectrum, is available only 
for non-polarized light and backscattering geometry on 
a polycristalline sample. To compare with experimental 
data, Eq. IKilmust be integrated over the solid angle by 
summing over all possible polarization vectors es and gl 
consistent with the backscattering geometry. The total 
cross section for unpolarized light in backscattering ge- 
ometry is obtained from Eq. 1161 with the substitution 

4 ( RI,x + Rly + RIz) + T{Rly + RIz + RIz) + 

( RxxRyy + RxxRzz + RzzRyy ) ^ l^S ' ' P 

(18) 

The resulting theoretical Raman spectrum is compared 
with experimental data in Fig. \7\ 



where n{uj) is the Bose factor, cos is the frequency of 
the scattered light, es and gl are the polarization vectors 
of the scattered and incident light, respectively. The Ra- 
man tensor R-' associated with the j-th phonon is given 
by 




(17) 



where Vo is the unit cell volume (274.25 A'^)^, r(K) is 
the position of atom K-th and x°° = (e°° — 5)/4tt is the 
electronic susceptibility. The inner longitudinal macro- 
scopic electric field has no effect on g-modes. The tensor 
RP is computed from by finite differences, by moving 
the atoms along the phonon displacement pattern with 
maximum displacement of 0.002 A. The Raman tensor 
(Eq. Il()|l for the Raman-active irreducible representations 
has the form 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 


E.<5) 

E.(3).E.(4),A (2) 

v.: J 


A,. (4) 
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FIG. 7: (Color online) Theoretical (black line) and experi- 
mental (gray line (red online), adapted from Ref.^) Raman 
spectra of polycristalline Mg(AlH4)2 for unpolarized light in 
backscattering geometry. The experimental spectrum is nor- 
malized by setting equal the experimental and theoretical in- 
tensity of the highest peak around 1800 cm~^. 



A 



Is 



The agreement is good, but for the highest frequency 
peak which is too high in energy and in intensity with 
respect to experiments. The displacement pattern of the 
Raman active intra-tetrahedra modes are sketched in Fig. 
0] The strongest Raman peak is due to the Aig stretching 
mode of the AI-H2 bond at 1845 cm~^. 



IV. CONCLUSIONS 



Ey,2 



— c —d 



— c 
-d 



The coefficients a,b,c and d calculated from first prin- 
ciples as outlined above are given for each mode in Table 
I. 



Based on Density Functional Theory we have opti- 
mized the structure of Mg(AlH4)2. Due to the lack of 
van der Waals forces within the current approximations 
to the energy functional, the interlayer spacing between 
the neutral A1H4-Mg-A1H4 sheets stacked along the c axis 
is largely overestimated (6 %) in our calculations. Con- 
versely, by fixing the lattice parameters to the experi- 
mental ones, the optimization of the internal structure 
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provides a geometry in fair agreement with that inferred 
experimentally from x-ray powder diffraction data. The 
main discrepancy regards the position of the hydrogen 
atoms (which however can not be detectd accurately from 
x-ray powder diffraction) resulting in Al-H bond lengths 
which differ up to 0.1 A from the exprimental ones. As 
a consequence, the AIH4 tetrahedra are much more sym- 
metric in the theoretical geometry than in that proposed 
experimentally. However, the IR and Raman spectra cal- 
culated within density functional perturbation theory are 
in good agreement with the experimental spectra which 



supports the correctness of the crystal structure emerged 
from the ab-initio calculations. 
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phonons within the Debye approximation as given above. 
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